A Further Analysis of a Cosmological Model with Quintessence and Scalar Dark 

Matter. 



Tonatiuh Matos* and L. Arturo Urena-Lopez^ 
Departamento de Fisica, 
Centro de Investigation y de Estudios Avanzados del IPN, 
AP 14-740, 07000 Mexico D.F., MEXICO. 
(February 1, 2008) 

We present the complete solution to a 95% scalar field cosmological model in which the dark 
matter is modeled by a scalar field $ with the scalar potential V($) = V [cosh (A ^/K^&) — 1] 
and the dark energy is modeled by a scalar field endowed with the scalar potential V("i/) = 
V [sinh (a ^/n^iS>)] . This model has only two free parameters, A and the equation of state w*. 
With these potentials, the fine tuning and the cosmic coincidence problems are ameliorated for 
both dark matter and dark energy and the models agrees with astronomical observations. For the 
scalar dark matter, we clarify the meaning of a scalar Jeans lenght and then the model predicts 
a suppression of the Mass Power Spectrum for small scales having a wave number k > k m in,<s>, 
where k m in,<s> — 4.5/iMpc -1 for A ~ 20.28. This last fact could help to explain the dearth of dwarf 
galaxies and the smoothness of galaxy core halos. From this, all parameters of the scalar dark 
matter potential are completely determined. The dark matter consists of an ultra-light particle, 
whose mass is m$ ~ 1.1 x 10" 23 eV and all the success of the standard cold dark matter model is 
recovered. This implies that a scalar field could also be a good candidate as the dark matter of the 
Universe. 

PACS numbers: 98.80.-k, 95.35. +d 



I. INTRODUCTION 

Observations of the luminosity-redshift relation of la 
Supernovae suggest that distant galaxies are moving 
slower than predicted by Hubble's law, implying an accel- 
erated expansion of the Universe jl| . These observations 
open the possibility to the existence of an energy com- 
ponent in the Universe with a negative equation of state, 
u> < 0, being p = tup, called dark energy. This compo- 
nente would be the currently dominant component in the 
Universe and its ratio relative to the whole energy would 
be J7a ~ 70%. The most simple model for this dark 
energy is a cosmological constant (A), in wich iv = — 1. 
The matter component VLm ~ 30% of the Universe de- 
composes itself in baryons, neutrinos, etc. and cold dark 
matter wich is responsible of the formation of struc- 
ture in the Universe. Observations indicate that stars 
and dust (baryons) represent something like 0.3% of the 
whole matter of the Universe. The new measurements 
of the neutrino mass indicate that neutrinos contribute 
with a same quantity like matter. In other words, say 
£Im = fl m + Qdm = Ob + fi„ + • • • + Qdm ~ 0.05 + Udm, 
where VLc d m represents the dark matter part of the mat- 
ter contributions which has a value of Qcdm ~ 0.25. 
The value of the amount of baryonic matter (~ 5%) is in 
concordance with the limits imposed by nucleosynthesis 
(see for example [||). Then, this model considers a fiat 



Universe (Qa + ~ 1) full with 95% of unknown mat- 
ter but wich is of great importance at cosmological level. 
Moreover, it seems to be the most succesful model fitting 
current cosmological observations || . 

However, from the theoretical point of view, A has 
some problems. First, the initial conditions have to be 
set precisely at one part in 10 120 , that is, an extremely 
fine tuning problem appears. Second, the cosmic coin- 
cidence: why is the current value of the energy density 
contribution of the cosmological constant of the same or- 
der than matter?. Third, particle theory predicts a zero 
cosmological constant, why is it not zero?. These prob- 
lems can be ameliorated by Quintessence, the model of a 
fluctuating, inhomogeneous scalar field (Q) rolling down 
a scalar potential V(Q) 0. It is assumed that flat mod- 
els with Hm — 0.33 ± 0.05 and a current value of the 
equation of state cuq = —0.65 ± 0.07 (uq can change 
along the evolution of the Universe) are the most con- 
sistent with all observations pj. However, there is not 
agreement about which scalar potential V(Q) is the cor- 
rect one. For instance, the pure exponential potential has 
been extensively analyzed [p|-|l"l| . It is known that there 
is a solution which makes the scalar energy density scales 
as the dominant background one, that could help to ame- 
riolate the fine tuning problem. Also, there is another 
solution that could make the Universe inflate, in good 
accord with SNIa observations. Moreover, in a scalar 
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dominated Universe, the scalar potential is effectively an 
exponential one A. But, nucleosynthesis constraints re- 
quire £Iq < 0.2, and then an exponential potential would 
never dominate the Universe |p|. 

Another example is a special group of scalar potentials, 
named tracker solutions M . The cosmology for these po- 
tentials is the same and independent of a large range of 
initial conditions (about 100 orders of magnitude), avoid- 
ing fine tuning. The equation of state uq changes with 
time towards —1 |^^|, then it can dominate the evolution 
of the Universe at late times. A typical example is the 
ure inverse power-law potential, V(Q) ~ Q~ a (a > 0) 
jig]. But the predicted value for the current equa- 
tion of state of the quintessence cannot be put in good 
accord with supernovae results Q. The same problem 
arises with other inverse power-law-like potentials. An- 
other possibility are the potentials proposed in |D| . They 
can solve the troubles stated above, but it is difficult to 
determine uniquely their free parameters. 

On the other side, we do not know the nature of 
the dark matter component Vic dm- The ACDM model 
agrees with the observations of large scale structure of the 
Universe, but cold dark matter over predicts subgalac- 
tic structure and singular cores for the halos of galaxies 
|l4| pUf . A scalar field model for dark matter has cre- 
ated a great expectation for solving the problem of the 
nature of dark matter too 13 17] pi] , A scalar field not 
only could give the correct energy density for the required 
matter in galaxies to predict the rotation curves of stars, 
but it is obtained the correct distribution of dark mat- 
ter in galaxies as well. Intriguinly, a natural solution is 
an exponential scalar potential jF7| . At the cosmological 
level, all attention has been put on the quadratic poten- 
tial <I> 2 , because of the well known fact that it behaves 
as pressureless matter due to its oscillations around the 
minimum of the potential ]22[| , implying that u>® ~ 0, for 
< p^ > — lu$, < p§ >, just like cold dark matter. 
In a recent paper p3] , we showed that the potential 
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is a reliable model for the dark energy, because of its 
asymptotic behaviors. At the other ha nd, a g ood model 
for dark matter could be the potential fll3| , pi|| 
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The mass of the scalar field $ is defined as m\ = 
U"|$ =0 = \ 2 k V . In this case, we deal with a massive 
scalar field. 

In this paper, we give explicitly all the solutions to the 
model at the cosmological scale and focus our attention 



on the scalar dark matter. In section II we find the com- 
plete solutions for a cosmology where the scalar field $ is 
the dark matter and the scalar field is the dark energy. 
We divide the evolution of the Universe in four stages: 
the radiation dominated era (RD) but before the oscil- 
lations of the scalar field 4>, the radiation dominated era 
with $ already acting as standard cold dark matter, the 
matter dominated era (MD) and finally the dark energy 
dominated era. We will recover the standard cosmolog- 
ical evolution of cold dark matter with the quintessence 
potential (|l|). In section III we will analyze the growing 
fluctuations of the scalar dark matter within the linear 
regime for the four stages of the cosmological evolution. 
From the equation of scalar fluctuations, we will deter- 
mine a Jeans lenght for scalar dark matter and show that 
modes larger than this Jeans lenght are growing modes. 
This last fact implies a cut-off in the Mass Power Spec- 
trum and we can fix all free parameters of the potential 
(||) by setting the cut-off wave number at the desirable 
value suggested in the literature. The dark energy fluc- 
tuations are not analyzed in detail because it has been 
already done in the literature. Finally, in section IV we 
summarize the results and give some future features to 
be investigated. For completeness, we will use previous 
results already shown in papers 21 2^| when necessary. 



II. SCALAR FIELD SOLUTIONS 

Due to current observations of CMBR anisotropy by 
BOOMERANG and MAXIMA Q, we will consider a flat, 
homogenous and isotropic Universe. Thus we use the flat 
Friedmann-Robertson- Walker (FRW) metric 



ds 2 



-dt 2 + a 2 (t) [dr 2 + r 2 (d9 2 + sin 2 O)d0 2 )] , (3) 



where a is the scale factor (a = 1 today) and we have 
set c = 1. The components of the Universe are baryons, 
radiation, three species of light neutrinos, etc., and two 
minimally coupled and homogenous scalar fields $ and 
^, which represent the dark matter and the dark energy, 
respectively. The evolution equations for this Universe 
are 



H z = [ - 

a 



$ + 3#$ 



dV(^) 
dV(^) 



= 



= 



P + 3H(p + p) =0, 



(4) 
(5) 

(6) 
(7) 



being k = 8irG and p (p) is the energy density (pressure) 
of radiation, plus baryons, plus neutrinos, etc. The scalar 
energy densities (pressures) are p$ = ^c|> 2 + U(<&) (p<s> = 
i $2 _ V(^)) and = i* 2 + U(*) (p w = i* 2 - )). 
Here dots denote derivative with respect the cosmological 
time t. 
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The numerical values that will be used along the paper 
are: the Hubble parameter, H a = 100 h km s Mpc~ x 
(H = 3.3 x 10" 4 hMpc' 1 in units of c = 1) with 
h = 0.65 ± 0.1, the current radiation energy density 
O7 /i 2 = 2.480 x 10~ 5 , the current baryon energy density 
&aBh 2 = 0.019 ± 0.0024, the background temperature 
T = 2.7277 ± 0.002 K and the amount of primordial 
helium Y He = 0.246 ± 0.0014 ||. We will consider 3 
species of light neutrinos with £1 „ = |(4/ll) 4 / 3 f2 07 per 
species. 



A. Radiation Dominated Era (RD) 

We start the evolution of the Universe at the end of 
inflation, i.e. in the radiation dominated era (RD). The 
initial conditions are set such that (pi$, Pi^) < Pi 7 . 

Let us begin with the dark energy. For the potential 
(Q) an exact solution in the presence of nonrelativistic 
matter can be found Jt],^3| . The parameters of the po- 
tential are given by 

— 3w* 







2^/3(1 + 
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with f2 D ^ and Q cdm the current values of dark energy 
and dark matter, respectively; and —0.6 > u>^ > —0.9 
the range for the current equation of state. With these 
values for the parameters (a, (3 < 0) the solution for the 
dark energy ("J) becomes a tracker one, is only reached 
until a matter dominated epoch and the scalar field 'J 
would begin to dominate the expansion of the Universe 
after matter domination. Before this, at the radiation 
dominated epoch, the scalar energy density p^ is frozen, 
strongly subdominant and of the same order than today 
p3| . Then the dark energy contribution can be neglected 
during this epoch. 

Now we study the behavior of the dark matter. For 
the potential (^) we begin the evolution with large and 
negative values of <i>, when the potential behaves as an 
exponential one. It is found that the exponential poten- 
tial makes the scalar field $ mimic the dominant energy 
density, that is, p$ — £i^p, A - The ratio of p$ to the 
total energy density is 



pa. 
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This solution is self-adjusting and it helps to avoid the 
fine tuning problem of matter, too. Here appears one 
restriction due to nucleosynthesis S 
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A 2 



< 0.2 
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acting on the parameter A > y24. We shall call a* the 
scale factor at the time when |Ay / «vI>| ss 1, i.e., when 
the scalar potential (||) is leaving the exponential behav- 
ior and entering into the polynomial one. The time evo- 
lution for RD before the scalar oscillations begin is (from 
eq- (§) 



(11) 



where H Q is the current value of the Hubble parameter. 
This result is the standard with an extra contribution 
due to the scalar contribution p$ to radiation. 

Once the potential (||) reaches its polynomial behavior, 
$ oscillates so fast around the minimum of the potential 
that the Universe is able only to feel the average values 
of the energy density and pressure in a scalar oscillation. 
During the time of a scalar oscillation, the term of the 
cosmological expansion in eq. (||) can be neglected, then 
we can write 



$ = . 



(12) 



Following p2| , from the time average in a period larger 
than a scalar oscillation but smaller than the Hubble time 
of the quantity 



~dt 



= $ 2 



it follows that 



A2 ^ dV 
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Taking only the quartic and quadratic term, as a good 
approximation, the average energy density and pressure 
can be expanded as |18] 



< P* > 

< P<s> > 
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then, the equation of state reads 

1 < $ 4 > 

< u<s> >= 



3 t&- < $ 2 > 



< $4 > 



(15) 
(16) 

(17) 



Observe that at the beginning of scalar oscillations, when 
the quartic term is still the dominant one in the poten- 
tial, the scalar field $ behaves as radiation like the early 
exponential behavior indicates. Once the quadratic po- 
tential dominates, we find that < <I> 4 >= (3/2) < $ 2 > 2 
and < p® >= 2 < V >. The equation of state then 
changes to 
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< UJ<s> > = 



16Vn 



< p* >, 



(18) 



with < cl>$ > going down to zero and < p$ > scaling 
as non-relativistic matter j2^] . A detailed analysis of the 
solution of the polynomial behavior can be seen also in 
fllf . Now, we would like the scalar field $ to act as cold 
dark matter, in order to recover all the successful features 
of the standard model. For this to be possible, we will 
first derive a relation between the parameters (Vb, A). 

If the transition from radiation to matter occurs 
smoothly (as suggested by eq. |l7])) at a — a*, then 



Pi* (a* 



PoCDM (a* 



Adjusting numerically and using eq. 
eq. ([l9l), we find that 

a 
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in the l.h.s of 
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On the other hand, now we use < p$ >~ Pcdm m the 
r.h.s of eq. @. Thus, 



V Q < A 2 k $ 2 > I „• 



PoCDM (a*) 
Therefore, the required relation reads pl[ 
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(21) 
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Notice that Vb depends on both current amounts of dark 
matter and radiation (including light neutrinos) and that 
we can choose A to be the only free parameter of potential 
(0). If we take A > 5, then we find the limit values 



K V > 10 6 Mpc~ 2 
m $ >5x 10 3 Mpc- 1 w 3 x 10" 26 eV 

XX C = (k V )- 1/2 <lkpc 
< w$ > - 10" 14 a~ 3 -> 0, 



(23) 
(24) 
(25) 
(26) 



where Ac = mj 1 is the scalar Compton length. Since 
now, we can be sure that p$ = Pcdm and that we will 
recover the standard cold dark matter evolution. More- 
over, it turns out that the scalar field <E> is an ultra-light 
cold dark matter particle. 

The dominant components of the Universe just after 
scalar oscillations begin are matter and radiation just 
like in the standard model. Then, we can give the time 
evolution with radiation and matter only for a > a* , ob- 
taining 
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where 
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is the scale factor at the time of radiation-matter equiv- 
alence, a 7 M ~ 3.3 x 10~ 4 (z ~ 3000). The standard case 
is recovered if a* — > (A — > 00) and that a* < a^M for 
A > 5, that is, the transition occurs before radiation- 
matter equality. 



B. Matter Dominated Era (MD) 

During this time, the scalar field <E> continues oscillat- 
ing and behaving as nonrelativistic matter and there is 
a matter dominated era just like that of the standard 
model. A short after matter completely dominates the 
evolution of the Universe, the scalar field \P reaches its 
tracker solution, eqs. (|J) and |p3fl 



<x v^o^ — arccoth 
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(31) 

l^V" =[(j3-l) tanlT 2 (a ^~ Q ^) + l] a 2 (3 V. (32) 



Po*a 

a (3 tanh -1 (a ^/kZ^) V, 



and it begins to be an important component. ajf$ is the 
scale factor at the time of equivalence between matter 
and dark energy, 



CLMW — 



a, 



M 



a, 



-1/(3^*) 



(33) 



where £l M and f2 0> i< are the current amounts of matter 
and dark energy, respectively. We find that 0.62 (z w 
0.6) > a M v > 0.73 (z w 0.37) for -0.6 > > -0.9. 

A tracker solution is recognized because the function 
r = {VV") I (V') 2 > 1 and nearly constant over the pos- 
sible initial conditions Q. From eqs. ([IT 32), 



r = [1 - /T 1 + /T 1 tanh 2 {a 



(34) 



thus r > 2 if lo^ < —0.6 for the possible initial con- 
ditions < $ < Mp. Therefore, potential ([!]) has a 
tracker solution. The time at wich the scalar field ^ de- 
pends upon initial conditions, but the late time behavior 
is independent of initial conditions over almost a range 
of 100 orders of magnitude 0,^3| . 

We can integrate eq. (|J) and then the time evolution 
for matter and dark energy only (neglecting radiation) 
reads MM 



H n t = 



|(i+w*) 
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2Fl (i£,£±i. 

2 1 2' 4' 4 ' 



Qm<B 



(35) 



4 



where 2F1 is the hypergeometric function. By analytical 
continuation, it can be shown that [E6l 



1F1 (u, v, w; z) = (1 — z) U iFi ( u,w — v,w 



z- 1 



If a -C fljj * in eq. (Ba) , then 



H t 
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Using eqs. (g|), it follows that 
/ 1 3 

2^1 -=i,T + i;i 



= l + o)*, 



therefore 



,3/2 



3 var 



A/ 



(36) 



(37) 



(38) 



(39) 



A complete numerical solution for the energy densi- 
ties p's and the dimensionless density parameters O's are 
shown in fig. || until today. The results agree with the 
solutions found in this section. It can be seen that eq. 
( p2| ) makes the scalar field 4> behave quite similar to the 
standard cold dark matter model once the scalar oscilla- 
tions begin and the required contributions of dark matter 
and dark energy are the observed ones J|l],p3| . 
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Let us try to give a complete solution for the time in the 
presence of radiation, matter and dark energy. Taking 
into account the limit value (|39|), we can stick up the 
solutions (^) and (|3^). Then, the complete time evolu- 
tion for a > a* (including the \& dominated era) can be 
written as 



FIG. 1. Evolution of the energy densities pi vs a: p 7 (pink), 
Pcdm (green), pb (red) and p<a (blue); with oj* = —0.7. 
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C. Scalar Field * Dominated Era (*D) 

At this time, the scalar field is the dominant com- 
ponente of the Universe and the scalar potential ([!]) is 
effectively an exponential one [p|]25|l. The time evolution 
becomes 



Hnt 



a |(l + W4,) 



3(l+o;$) y/n^ 



(41) 



thus the scalar field if? drives the Universe into a power- 
law inflationary stage (a ~ t p , p > 1). Note that usual 
tracker solutions do not have this late exponential be- 
havior. Also, tracker equation of state usually changes 
toward —1 once the scalar field dominates finishing the 
latter as a cosmological constant Q . But the scalar equa- 
tion of state o;* will never change its tracker value and 
for potential (|l|) there is no solution with oj* = — 1 either 



FIG. 2. Evolution of the dimensionless density parameters 
vs the scale factor a with £1 m = 0.30: ACDM (black) and 
*$J)M for two values of A = 6 (red), A = 8 (green). The 
equation of state for the dark energy is o;* = —0.8. 

The current age of the Universe can be calculated from 
eq. © with o = l. Thus, 1.09 > H t > 4.02 for 
—0.6 > 0;* > —0.9. The age of the Universe could be 
(15 < t a < 60) x 10 9 years, irrespectively of the value of 
A. We recovered the standard cosmological evolution and 
then we can see that potentials (0^) are reliable models 
of dark energy and dark matter in the Universe. 



III. LINEAR PERTURBATION THEORY 

In this section, we analyze the perturbations of the 
space due to the presence of the scalar fields \& using 
analytical approximations. First, we consider a linear 
perturbation of the space given by hij. We will work in 
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the synchronous gauge formalism, where the line element 
is ds 2 = a 2 [— dr 2 + (Sij + hij)dx t dx J ]. The equations for 
perturbations in the £;-space (k — kk) are 

hij {x , t) = 



d 3 ke lk - s 



(fc,r) + [^kikj - -Sij )<;,/(/.-. t)\. (42) 



A. Scalar Perturbations 

We must add the perturbed equations for the scalar 
fields $(r) -> $(r) +^(fc,r) and *(r) -> *(r) +^(fc,r) 



5/3* = -^$0 + V'<f>, 



k r\ 
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ft = 47rGa 2 (5r n 
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ft + 2-ft - 2k 2 r] = -8ira 2 6Ti, 
a 



(43) 



k 2 f 1 = ^TTGa 2 (p + p)9, (44) 



(45) 



ft + 6^ + 2- (^ft + 6?7j -2k z r)= -2AnGa z {p + p)cr- (46) 

where ft is the trace of the metric perturbations hij ; 6T£ 
is the perturbation to the momentum-energy tensor and 
now dots denote derivatives with respect to conformal 
time r. The velocity 9 and shear a perturbations are 
defined as 



5p$ = ^r<j>0 - V'0j 



1 



(56) 



5p* = \vif) + V'tp 
Spy = Xr^tp - V'tp 

with the evolution equations for the perturbations 



(57) 
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kikj ~^&ij ) Sj'j 



(47) 
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The density contrast (8 = dp/ p) wich accounts for the 
over density relative to the homogeneus cosmological en- 
ergy density and 9 for perfect fluids are governed by 



1 



2H<p + k z (f> + a A V"<p + -$ft = 
1 • : 



(58) 



i> + 2Hip + k 2 ip + a 2 V"'>p+ = 0. (59) 

Here, primes are derivatives with respect of the unper- 
turbed scalar fields <& and 'J', respectively. Before we 
solve these equations, we analyze first the meaning of a 
scalar Jeans lenght. 



W^PS-**. (51) 



The evolution equations with the unperturbed FRW met- 
ric (ky-[7f) in terms of the conformal time are 



H 2 
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- 2H^> + a 2 
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p + 3H(p + p) 


= o, 


(55) 



B. Damping of the Scalar Power Spectrum 



Recalling the results for the average scalar pressure 
and energy density (eqs. 15 1^), the velocity of sound in 
the scalar fluid $ is 



1 

8Vo 



< >~ Kr U eT 3 



(60) 



the velocity of sound decreases rapidly to zero as w$ does. 
As suggested in the literature, we could define an effective 
Jeans length [OJllpl], 



G<p> 



= 7T {k V ) 



-1/2 



(61) 



Ti. being the conformal Hubble factor, the scalar en- 
ergy densities are p$ = (l/2a 2 )$ 2 + V($) and p^ = 
(l/2a 2 )^> 2 + V (/<&). The scalar pressures are p$ = 
(l/2a 2 )$ 2 - V($) and p* = (l/2a 2 )* 2 - V(V). The 
solution to this equations are those already found in the 
previous section when written as functions of the scale 
factor a 



from which we obtain that Lj < 3.14 kpc. Therefore, 
we can expect that the scalar field fluctuations produce 
only complex stelar objects bigger that this Jeans length. 
In first approximation, this model could explain the sup- 
pression of subgalactic structure p3[ . Also, it is likely 
that the quantity r c = (kqV^)" 1 / 2 could be important at 
galactic scale [gl|. But, in order to have a good estimate 
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of the real Jeans lenght, we need to analyze the flutu- 
ation equations. We will find that eq. ([3l|) is only an 
approximation to the real thing. 

It is known that scalar perturbations can only grow if 
the &; 2 -tcrm in cqs. ( |5^j5^ ) is subdominant with respect 
to the second derivative of the scalar potential, that is, if 
k < aV" |2^] . Let us anlyze the scalar dark matter case. 

According to the solution given above for potential (§) , 
the behavior for = aV" is 



/,-„(„) = < 2H a ^^Xa 1 (a<a*) (g2) 

(a* < a) 



with a minimum value given by 



= m$a* ~ 1.3 Ay A 2 — 4 



&0CDM 



H , (63) 



and then for A > 5, k m i n ^ > 0.375 Mpc' 1 (see fig §. 
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FIG. 3. Evolution of the Hubble wave number ku (red) 
and the scalar wave numbers fc* (green) and fc* (blue). 

It can be assured that there are no scalar perturbations 
for k > ra$, that is, bigger than k$ today. These k corre- 
sponds to scales smaller than 1.2 kpc (here A = 5). They 
must have been completely erased. Besides, modes which 
m$ > k > femin,* must have been damped during certain 
periods of time determined by 



2iJ A 



O k 

X 2 - 4 m$ 



(64) 



From this, we conclude that the scalar power spectrum 
of $ will be damped for k > k m i n ^ with respect to the 
standard case. Also, in fig. || we can see that the Hubble 
wave number ku is greater than fc$ until a little after a* . 
It means the Jeans length (|65| ) is well within the Hubble 
horizon for most of the time, then the scalar fluctuati- 
nos can form structure at scales smaller than the Hubble 
lenght H . 

Therefore, the real Jeans length must be 



Lj{a) = 2irk„L 9 , 



(65) 



and it is a universal constant because it is completely 
determined by the mass of the scalar field particle. Lj 



is not only proportional to the quantity (kqVq)~ 1 / 2 , as 
suggested by the approximation in eq. (pl|), but also the 
time when scalar oscillations start (represented by a*) is 
important, as it was suspected in jl3|. 

On the other hand, the wave number kq, = aV" is al- 
ways out of the Hubble horizon, then only structure at 
larger scales than H -1 can be formed by the scalar fluc- 
tuations Instead of a minimum, there is a maximum 
k m ax^! ~ 10 -3 . All scalar perturbations of the dark en- 
ergy which k > 10~ 3 must have been completely erased. 
Perturbations with k < 10~ 3 have started to grow only 
recently (see fig. ||). For a more detailed analysis of the 
dark energy fluctuations, see pg,B(|. 



C. Scalar Power Spectrum for $ 

In this subsection, we will find some analytical approx- 
imated solutions to eq. (j58|), having in mind the Jeans 
Length (|65|), that is, we will analyze eq. ( |5S| ) only for 
modes k < k m in,<s>- 

We can have some physical insight into the previous 
differential equations if we write the evolution equations 
(E1EBI) in the form 



= K a 3 {ST°-5Tj 



(66) 



then we can change this into 



d_ 

5p$ 



(ah) 



5p<s> 



OA 1 + + fM* 1 + 3-^- 



5. 

Spy 



b°b- 



= 0, (67) 



where Oj = pi/pr, Pt = {SH 2 ) / (K n a 2 ) . For perfect flu- 



= v s ~ uj. For scalar fields, 
in general and so we have 



ids, it happens that (Sp/Sp) 
we can not identify v 2 
explicitly written (<5p$^ /5p<s>^)- 

We can see easily that the dominant background com- 
ponent of the Universe is the dominant term in the differ- 
ential equation (^), too. Thus, we can use the standard 
results for the fluctuations during ceratin stages of the 
cosmological evolution. Also, it is worth to mention that 
all perturbed quantities can be written in terms of the 
trace of the metric perturbations h, then it is enough to 
have a solution for it. Sec (27],|3l[] for the other equations 
required and how to solve them. 



1. a < a* : Radiation Dominated Era 

During RD, the scalar energy density /o$ evolve as a 
perfect fluid with constant equations of state o;$ = 1/3 
(remember that the exponential potential mimics the 
dominant energy). For the scalar field $ we have 
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= 

a 2 V' = 



4 4 



4 

A K 



(68) 



(69) 



a 2 V" 



AH 2 



Since ri = r 1 and radiation dominates equation (67), 
/i evolves as in the standard case, h — Ct for modes 
out the Hubble horizon and h — Ct" 1 for modes inside 
the Hubble horizon |}2| (for a detailed calculus with the 
exponential potential, see Q), being C = const.. 

Having that a = r, the evolution equation rt58) for a 
mode out of the Hubble horizon can be written (see eqs. 



2 

-< 

T 



2C 



(70) 



Thus, the growing solution for <p is 

C , 



k o 0(t) 



5A 



(71) 



where we can recover the result — (l/2)/i = 8c dm = 
(5A/4) v /k^> ||. Here, 8c dm would be the standard 
CDM density contrast. This leads to 



<5* — -rpdcDM, 
15 



(72) 



that is, the scalar density contrast 5$ evolves as the stan- 
dard one but with smaller amplitude. Note that this re- 
sult is independent of A. This result can also be obtained 
from eq. (pQ) with 



8p<s> 



3, 



I 

3' 



(73) 



where we can observe that the scalar field <!> does not 
behave completely as perfect fluid. 

For modes inside the Hubble horizon, the equation to 
be solved is 



T 



2C 1 

Kq A T 2 



(74) 



The general solution is of the form 



C 



k 4>(t) = -— = const. 



(75) 



and then 8$ = (C/6) = const. The modes inside the 
Hubble horizon do not grow during RD. 



2. a* < a: Radiation and Matter Dominated eras 



Once the scalar field <& begins to oscillate, it happens 



that V" = 



Recalling that k < k 



s) can be written 



1 



+ 2H(j) + CT7740 + 



o. 



< a 2 V" 



eq. 



(76) 



The scalar perturbation <f> oscillates with the same fre- 
cuency that the unperturbed Following the previous 
section, if we take the time average of the quantity 



d($(j>) 
dr 



(77) 



we obtain 



< §<j) >= -<$</>> 



(78) 

The second and fourth terms of eq. (|7^) are almost con- 
stant during the time of a scalar oscillation. Then, 



< $</> >« - < a 2 m|,< 
Therefore, we find that 

< <5p$ > <<£(/>> 



>= - < a 2 V'c 



< a 2 V'<p > 



> 



(79) 



(80) 



< Sp® > <$(/>> + < a?V'4> > 

It is now convenient to rewrite eq. (|7^) in the same form 
as eq. (|50|), 



< <5p$ > < p$ > 

< Sp® > < p<s> > 



= - 1 



< pg> > A h 



< pg, > J 2 

(81) 

Therefore, with < w$ > going to zero (eq. |l8|) it follows 
<S$ = S C dm- (82) 

Due to its oscillations around the minimum, the scalar 
field $ changes to a complete standard CDM and so do 
its perturbations. All the standard growing behavior for 
modes k < fcmin,* is recovered and preserved until to- 
day by potential (||). In fig. [|, a numerical evolution of 
5 = Sp/p is shown compared with the standard CDM 
case [^IJ. The results ( |72] , |82| ) agree with the numerical 
solution. The numerical evolution for the density con- 
trasts was done using an amended version of CMBFAST 

M. 
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FIG. 4. Evolution of the density contrasts for baryons 8t, 
standard cold dark matter 8c dm and scalar dark matter 8<s> 
vs the scale factor a taking S1 m = 0.30 for the models given 
in fig. (^|). The modes shown are k = O.lMpc -1 (top) and 
k = 1.0 X 10~ 5 Mpc' 1 (bottom). 

In fig. ^| we can see S 2 at a redshift z = 50 from a com- 
plete numerical evolution using the amended version of 
CMBFAST. We also observe a sharp cut-off in the pro- 
cessed power spectrum at small scales when compared 
to the standard case, as it was argued above. This sup- 
pression could explain the smooth cores of dark halos in 
galaxies and a less number of dwarf galaxies JTa] . 




k (h/Mpc) 

FIG. 5. Power spectrum at a redshift z = 50: ACDM 



(black), and $CDM with A = 5 (red) and A 
The normalization is arbitrary. 



10 (blue) 



At this point, we would like to mention some co- 
incidences between some results in ref. |2p| l in wich a 
quadratic potential is used and the results of this sec- 
tion. In that reference, it is argued that the Jeans lenght 
is the de Broglie wavelenght at the ground state of the 
particle in the gravitational potential well and that the 
power spectrum is supressed relative to the CDM case. 
We see before that most of the interesting properties of 
potential (0) as dark matter are due to its polynomial 
behavior $r. Then, it is not strange that, in our case, 
the mass power spectrum is also related to the CDM case 
by (see M) 



P*(fc) 



1 + x b 



P, 



CDM 



(k), 



(83) 



but using x — (fc/fcmin,*) with k m i n ^ being the wave 
number associated to the Jeans lenght (65). The differ- 
ence with respecto to the case treated in |20(] is that the 
relevant time scale is that when scalar oscillations start 
and not that of radiation-matter equality. 

If we take a cut-off of the mass power spectrum at 
k = 4.5/iMpc -1 pi , we can fix the value of parameter 
A. Using eq. (|63|), we find that 

A ~ 20.28, 

V ~ (3.0 x 10~ 27 Mpi ~ 36.5 eV) 4 , (84) 
1.1 x 10^ 23 eV. 



to* ~ 9.1 x 10" 52 Mpi 



where M Pl = 1.22 x 10 19 GeV is the Plank mass. All pa- 
rameters of potential (||) are now completely determined 
and we have the right cut-off in the mass power spectrum. 



3. Scalar Field ^ Dominated Era 

Only for completeness, we will draw some general fea- 
tures of the evolution of fluctuations during the dark en- 
ergy dominated era. 

At this era, the scalar field \P now dominates both the 
evolution of the Universe and the differential equation 
(§7|). We do not worry anymore of $, because its per- 
turbed solution continues being <5$ = —{l/2)h due to its 
oscillations around the minimum of the potential. The 
scalar energy p$ evolve as a perfect fluid with equation 
of state tuqt due to the effective exponential behavior of 
potential (|l|). Then, 



V 



a/1 + \J a 2 p^ 

1 — U)^ 
7, 



(85) 
(86) 
(87) 



V' = y/Z{l + UJy)V 

V" = 3{1 + uj^)V. 



Since the scalar field \& dominates the evolution of the 
Universe, it is straightforward that 



m 2 

K 



n = 



1 + 3w$ 



(89) 
(90) 



being Too the size of the event horizon (r — > t^) p7[ . 
Since we are interested in possible growing modes, we 
will consider eq. (|59|) only in the case k 2 <C a 2 V" . The 
evolution equations ([39], [37|) become: 



tp + 2B 
h + B 



c- 



D 



(too - r) " (too - t) 2 2 (too - t) 
h . _ ip %j) 



{Too ~ T) 



AD 



(Too - T) E '(t qo - T) 



= (91) 
2 = (92) 
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where the constant coefficients are 
2 



B 



1 + 3w$ 



c = |(i-c4)s 2 



D = ^3(1 + 



£ = -V3(1 + ^*)(1-^*)B 2 - 



(93) 

(94) 
(95) 
(96) 

-1/2 



The scalar fields were normalized in units of n t 
We can try solutions of the form h = h (r — T 00 ) m , 
tp = tp (r — Too)™, where the values of m are the so- 
lutions to the equation 

Dtyi 

m [m(m - 1) + 2Bm + C] [(m - 1) +B] = — {ADm - E) 

(97) 

with an obvious solution m = and two complex roots 
. The other real root gives m > if —0.6 < w<p < 
-0.9. For instance, if u^, = -0.6, m = {7.83, 0.835 + 
3.2i, 0.835 — 3.2i}. Then, all of the solutions are decay- 
ing ones. This result is not surprising because by this 
time the Universe has already entered in an inflationary 
stage. Then, we conclude that perturbations in the linear 
regime will be suppressed by the accelerated expansion 
of the Universe. 



IV. CONCLUSIONS 

We have developed most of the interesting features of 
a 95% scalar-nature cosmological model. The interesting 
implications of such a model are direct consequences of 
the scalar potentials (j^). 

On one hand, we have modeled the dark energy of 
the Universe using a scalar field with a sinh potential 
|p3| . This potential has the advantages that at early 
times it is a pure inverse power-law one, V{^) ~ , 5~ Q! , 
thus the cosmology at late times is extremely insensi- 
tive to initial conditions, reducing the fine tuning. The 
'J solution is only reached during a matter dominated 
phase of the Universe, thus it is natural the existence of 
a current dark energy dominated epoch; the cosmic co- 
incidence is ameliorated. Its late exponential behavior 
drives the Universe into a power-law inflationary stage, 
in good accord with supernovae results. Nevertheless, 
a fine tunning problem remains in determining the pa- 
rameters of the potential. Also, we do not know about 
a fundamental theory that could predict this kind of po- 
tential. However, quintessence models with a expectation 
value of the field of the order of the Planck mass can be 
considered within supergravity |p3[| . 

On the other hand, we have modeled the cosmologi- 
cal dark matter using another scalar field $ with a cosh 
scalar potential. As we have shown in this work, the solu- 
tions found alleviate the fine tuning problem for cold dark 



matter, too. Once the scalar field $ begins to oscillate 
around the minimum of its potential (J2J) , we can recover 
the evolution of a standard cold dark matter model be- 
cause the dark matter density contrast is also recovered 
in the required amount. It should be noticed that the 
results (about growing density perturbations) are inde- 
pendent of the parameters of potential (Q) . Thus, the pre- 
dicted angular and power spectrums in the linear regime 
are those already shown in p3| . These spectrums are 
subject only to the imprint of the scalar potential (Q). 

We also find an important difference with respect to 
the standard dark matter model. Analyzing the fluc- 
tuation equations, we clarified the meaning of a Jeans 
length for this model: it is related to the mass of scalar 
particle and to the time when scalar oscillations start. 
This Jeans lenght provokes the suppression in the power 
spectrum for small scales, that could explain the smooth 
core density of galaxies and the dearth of dwarf galaxies. 
Up to this point, the model has only one free parameter, 
A. However, if we suppose that the scale of suppression 
is k = 4.5/iMpc, then A w 20.28, and then all param- 
eters are completely fixed. From this, we found that 
Vo ~ (36.5 eV) 4 and the mass of the ultra-light scalar 
particle is to$ ~ 1.1 x 10~ 23 eV. The quantity (k Vo) -1 / 2 
or the scalar mass m$ could play an important role in 
galaxies, appearing maybe in the observed constant core 
density of dark halos (see [ pT)Jl^ ] ) . This last fact could be 
a signature of the parameters of the dark matter potential 
(Q) as it has been shown that an exponential potential ap- 
pears when analyzing scalar dark matter at galactic level 
p7| . Further investigation will be published elsewhere. 

A question could arise here: How good is eq. (^2|)?. 
Why does f2o 7 also appear?. A possible answer could be 
that, in fact, potential (|^) can be written 



U [cosh(A v ^$)-l] = 2U 



(98) 



Then, potential (^) is another sinh-like potential. Ob- 
serve the similarity between eq. ( |22| ) and the last of eqs. 
(||), being the last one a generic solution of sinh-like po- 
tentials W: They both involved the previous dominant 
component. And it is this singular feature that make us 
think about a cosmic coincidence of matter. Therefore, 
it would be clear that after a radiation-dominated era, 
a matter-dominated era must appear. Moreover, after 
a matter-dominated era, a dark energy-dominated era 
must appear, too. 

Summarizing, a model for the Universe where the dark 
matter and energy are of scalar nature can be realistic 
and could explain most of the observed structures and 
features of the Universe. 
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